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Abstract 

In this paper we propose a Hamiltonian generalizing the interaction of the two- 
level atom and both the single radiation mode and external field • • • a kind of cavity 
QED. We solve the Schrodinger equation in the strong coupling regime by mak- 
ing use of rotating wave approximation under new resonance conditions containing 
the Bessel functions and etc, and obtain unitary transformations of four types cor- 
responding to Rabi oscillations which perform quantum logic gates in Quantum 
Computation. 

In this paper we consider a unified model of the interaction of the two-level atom 

and both the single radiation mode and external field (periodic usually) in a cavity. We 

deal with the external field as a classical one. As a general introduction to this topic in 

Quantum Optics see pQ, |2j, [H|. Our model is deeply related to the one of trapped ions 

in a cavity with the photon interaction (Cavity QED). 
*E-mail address : fujii@yokohama-cu.ac.jp 



In our model we are especially interested in the strong coupling regime, [5], [E], [T5] . 
One of motivations is a recent very interesting experiment, [T2~] . See [3] and JU] as a 
general introduction. 

In [0] and [Hj we have treated the strong coupling regime of the interaction model of 
the two-level atom and the single radiation mode, and have given some explicit solutions 
under the resonance conditions and rotating wave approximations. 

On the other hand we want to add some external field (like Laser one) to the above 
model which will make the model more realistic (for example in Quantum Computation). 
Therefore we propose the unified model. 

We would like to solve our model in the strong coupling regime. Especially we want to 
show the existence of Rabi oscillations in this regime because the real purpose of a series 
of study ([Sj, [Hi], [TS] ) is an application to Quantum Computation (see [T3] as a brief 
introduction to it). 

In this paper we solve the Schrodinger equations in this regime by making use of rotat- 
ing wave approximation under new resonance conditions containing the Bessel functions 
and etc, and obtain a unitary transformation in one qubit case and unitary transforma- 
tions of four types in two qubit case, which perform quantum logic gates. 

Our solutions might give a new insight into Quantum Optics or Condensed Matter 
Physics as well as Quantum Computation. 

Let {o"i, 02, 03} be Pauli matrices and I2 a unit matrix : 
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and a + = (1/2) (oj. + i&i), o\_ = (1/2) (01 — i<T\). Let W be the Walsh-Hadamard matrix 




(2) 

then we can diagonalize o\ as o\ = \Va3W~ 1 = Wa^W by making use of this W. The 
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eigenvalues of o\ is {1, — 1} with eigenvectors 



V 





(3) 



1 
1 

Let us consider an atom with 2 energy levels Eq and E\ (of course E\ > Eq). Its 
Hamiltonian is in the diagonal form given as 

I E 

Ho=\ |. (I) 

V Ex 

This is rewritten as 

/ 1 n \ „ ~ / 



_E + -^1 
-"0 = 7, 



V 




1 \ Ek-Eo 1 A 

= A 1 2 - — cr 3 , (5) 



1 

where A = E% — Eq is a energy difference. Since we usually take no interest in constant 
terms, we set 

#o = ~<r 3 . (6) 

We consider an atom with two energy levels which interacts with external (periodic) 
field with gcos(u>Et)- In the following we set H — 1 for simplicity. The Hamiltonian in the 
dipole approximation is given by 

A 

H = H + gcos(u E t + <p)a! = - —a 3 + g cos(u E t + (pja^ (7) 

where uje is the frequency of the external field, g the coupling constant between the 
external field and the atom. We note that to solve this model without assuming the 
rotating wave approximation is not easy, see |Hj, [H], |2I], |22] |2B]- 
In the following we change the sign in the kinetic term, namely from —A/2 to A/2, to 
set the model for other models. However this is minor. 

Now we make a short review of the harmonic oscillator within our necessity. Let a{a*) 
be the annihilation (creation) operator of the harmonic oscillator. If we set iV = a) a (: 
number operator), then we have 

[N, a 1 "] = a f , [N, a] = -a , [a f , a] = -1 . (8) 



Let H be a Fock space generated by a and a\ and {\n)\ n e N U {0}} be its basis. The 
actions of a and a' on 7i are given by 

o|n) = \/n|n — 1) j a |w) = \/n + l|n + 1) , iV|n) = n|n) (9) 

where |0) is a normalized vacuum (a\0) = and (0 1 0) = 1). From (JUJ) state |n) for n > 1 
are given by 



\n) 



10) 



(10) 



These states satisfy the orthogonality and completeness conditions 

oo 

(m\n) = 5 mn , ^2\n)(n\ = l. (11) 

n=0 

Then the displacement (coherent) operator and coherent state are defined as 

D(z) = e za] ~~ za ■ \z) = D(z)\0) for z e C. (12) 

We consider the quantum theory of the interaction between an atom with two-energy 
levels and single radiation mode (a harmonic oscillator). The Hamiltonian in this case is 



H = co>l 2 <8> a) a + — er 3 (8)1 + go\ ® (a^ + a) 



(13) 



where u is the frequency of the radiation mode, g the coupling between the radiation field 
and the atom, see for example jH], jH]- 

Now it is very natural for us to include (JZJ) into (JESJ), so we present the following 



General Hamiltonian 



A 



H = uj! 2 ® a) a + g\G\ ® {a} + a) + — cr 3 (g) 1 + g 2 cos(uEt + 0)(Ti <8> 1. (14) 

Our Hamiltonian has two coupling constants. We note that our model is deeply related 
to the Cavity QED (trapped ions in a cavity with the photon interaction). 

10} 
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This Hamiltonian is also related to the one presented recently by Schon and Cirac 0] 



H 



2m 



+ lo 1 2 <g> a) a + g(x) (a + (g> a + a__ <g> a)^j + 

y <x 3 ® 1 + ^ (e"^' (T+ g> 1 + e iWii <7_ ® l) . 



(15) 



For the meaning of several constants see [1]. They have assumed the rotating wave 
approximation (see for example [3]) and the resonance condition, and use a position- 
dependent coupling constant g(x), so their model is different from ours in these points. 
A comment is in order. Following |lj the Hamiltonian dHj) might be modified to 

p 2 A 

H = h ul 2 <8> a? a + gi(x)a 1 ® (a* + a) + —0-3 <g> 1 + g 2 cos{uj E t + 0)<Ti <g> 1. (16) 

2m 2 

This model is a full generalization of (|15|). however we don't consider this situation in the 

paper. 

We have one question : Is the Hamiltonian ()14j) realistic or meaningful ? The answer 
is of course yes. Let us show one example. We consider the (effective) Hamiltonian pre- 
sented by NIST group [3] , jHj which were used to construct the controlled NOT operation 
(see [T3| as an introduction). 



H = W0I2 <S> a) a + g C 



a + ® e iT?(at+a) + a_ 



-ir](a^ +a) 



+ — <x 3 <g> 1. 



(17) 



We can show that under some unitary transformation the Hamiltonian (|17jl can be 
transformed to f)14j) with special coupling constants, [7j. This is important, so we review 
and modify 

We set 2 A = irj(a> + a) for simplicity, then 




=U(rj)(<r 3 ®l)U(riy 



18) 



where 

/ 



1 / e A \ 1 1 \ 1 



e A — e A 



e A e" A 



(19) 



and e A = D(ir]/2) where D(0) is a displacement (coherent) operator defined by p2j) . 
Then it is not difficult to show 

U{r]) ] HU{r]) = ^^1 2 ®1 + uj q 1 2 ® a) a+^-ax® {-ia ] + ia)+ga^®l-^a 1 ®l. (20) 

To remove i in the term containing a we moreover operate the unitary one 



(l 2 ® e l{n/2)N ) U{rtfHU{r)) (l 2 ® e< v /*> N 



= ^^1 2 ®1 + ® « T a + -^o-i ® K + a) + ® 1 - — ai <g> 1, (21) 

~t ^ — 

where we have used the well-known formula 

e l0N ae~ i0N = e~ l9 a, e^aV*" = e l V, 

see [TE]. Since U{rj) can be written as Uirf) = (cr + <g> e A + cr_ <g> e~ A ) (W ® 1) , so if we 
write 

T{rj) = C/(t7)(1 2 ® e"^ /2)Af ) = (<x + ® e A + a_ ® e - A )(W^ ® e" i(7r/2)7V ), (22) 
then we have 

T{r])^HT{ri) = ^p-l 2 ® 1 + ® a f a + <g> (a f + a) + g<r 3 <g> 1 - ^ <g> 1. 

4r ^ 

Here we have no interest in the constant term, so we finally obtain 

H = T( V ) |cj 1 2 ® a ] a + ^-a t ® (a f + a) + ga 3 ® 1 - —a x ® l| T(^) t . (23) 

T(?y) is just the unitary transformation required. We note that the last term is constant, 
which case is a special one. 

At this stage we would like to make a further generalization of the Hamiltonian (fUjl 
to make wide applications to Quantum Computation.® 
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Let {K + , K_, K 3 } and {J + , J_, J 3 } be a set of generators of unitary representations 
of Lie algebras su(l, 1) and su(2). They are usually constructed by making use of two 
harmonic oscillators (two-photons) a x ,02 as 

su(l, 1) : K + = a^a 2 \ K_ = a 2 a x , K 3 — - (a^ai + 0,2*0.2 + l) , 
su{2) : J + = 0^02, J- = a 2 ^a x , J 3 = - (a x Jf a x — 0,2^0-2) • 

Then we can make the similar arguments done for the Heisenberg algebra {a\a,N}, 
namely (JHJ ~ (fT^jl . see for example [TH] . 

We have considered the following three Hamiltonians in 



(N) H N = ul 2 ®a)a + ^a 3 ®l + ga x ®{a ] + a), (24) 

(K) H K = ujl 2 ®K 3 + ^-a 3 ® 1 K + ga x ®(K+ + K_), (25) 

(J) Hj = lu1 2 ® J 3 + y^s ® 1j + ^1 ® (J+ + J_). (26) 

To deal with these three cases at the same time we set 



(N) {a\a,N}, 

(K) {K + ,K.,K 3 }, (27) 
(J) {J + ,J-,J 3 } 



and 



H L = uj1 2 ® L 3 + y(T 3 ® 1 L + 0<ti ® (L+ + L_). (28) 



Therefore the Hamiltonian that we are looking for is 
Unified Hamiltonian 

H L = ul 2 ®L 3 + g x a x <g> (L+ + L_) + y<x 3 ® 1 L + g 2 cos(u E t + (j))a x <g> 1 £ . (29) 

By the way, from the lesson in ()23|) we can also consider a weak version of (j2Uj) with 
the last term being constant 

H L = uj1 2 <g> L 3 + ^i(Ti <g> (L+ + L_) + ^-cr 3 <g> 1 L + g 2 a 1 ® 1l, 



where we have reset g 2 = g2cos(0) for simplicity. This is not a bad model as seen from 
(Unj). This restricted model has been treated in |T8] . 

Now we would like to solve the Hamiltonian (J2HI), especially in the strong coupling regime 
(9i > A). 



Let us transform (ffijj) into 

H L = 1 2 ® ujL 3 + <7i® {pi(L + + L_ 
A 

= H + —cr 3 <g> 1 L . 



The method to solve is almost identical to 
use of the Walsh-Hadamard matrix (12*1) 



g 2 cos(u E t + 0)1l} + ycr 3 <8> 1l 



(30) 



so we give only an outline. By making 



(Q, x) 



31) 



F = (W ® 1 L ) [1 2 <g> wL 3 + 3 ® {oi(L+ + L_) + g 2 cos(cu E t + <j>)l L }] (W 1 <g> 1 L ) 
= E (|A)®e-^ (L+ - i - ) ){fiL 3 + Ap 2 cosKt + 0)l L }((A|®e^ (L+ - L - ) 

A=±l 

where |A) is the eigenvectors of 0\ defined in Q and Q, x are given as 

(N) u, x = 2g 1 /u, 

(K) uj^I - (2 9l /ujy, x = tanh- 1 (2p 1 /c^), 
[J) uyjl + [2 gi /uj)\ x = tan- 1 (2p 1 /w). 

That is, we could diagonalize the Hamiltonian Hq. Its eigenvalues {E n (t)} and eigenvec- 
tors {|{A,n})} are given respectively 

(iV) n(-&+n) + \g 2 coa(w E t + </>), |A) O e"^^-^ 

(K) Q(K + n) + Xg 2 cos(u E t + (/>), \\) ® e-^ K +- K ^\K,n), 

(J) 0(-J + n) + Xg 2 cos(u E t + 0), |A) ® e ~^ (J +~ J - ) | J,n) 

(32) 

for A = ±1 and n6NU {0}, where E n (t) = E n + Ap2Cos(o;£;t + 
be written as 

Ho = T,T, E n(t)\{\n})({X,n}\. 



(E n (t), |{A,n})) 



Then Hq above can 



A n 



Next we would like to solve the following Schrodinger equation 



(33) 



To solve this equation we appeal to the method of constant variation. First let us solve 
i4;ty = H ^/, which general solution is given by = Uo(t)^Q, where \l/o is a constant 
state and 

u (t) = Y,T, e ~ i{tEn+x{92/WE)sm{ " Et+ ^ )} \{^ n })({^ n }\- ( 34 ) 

A n 

The method of constant variation goes as follows. Changing like \l/o — > *&o{t), we have 

ijJS>Q = y^oVs ® WUo$o = y^F^o (35) 
after some algebra. We must solve this equation. H F is 

H F = ^^ e ft(Bra -^ )+<(A -' i)(fla/c,,B) ' <n(WBt+ * ) ({A,m}|(t73 ® 1l)|{m,™}> |{A,m})({//,n}| 

= ^^ e i{(%-n) + 2Afc/, E ) sm (^)}^^| e Ax(L + -L_)| n ^ | {A, 777,}) ({- A, 7i} | , (36) 
A m,n 

where we have used (A|o"3 = (— A| and |n)) is respectively 



In the following we set for simplicity 



(N) |n), 
(K) \K,n), 
(J) |J,n). 



9(0 = , 2 sin( ^ t + ^ . (37) 
we 

Here we divide H F into two parts H F = H F + H F where 

H'f = EE^({n|^- L )|n)) |{A,n})({-A,n}|, (38) 

A n 

#f = E E e i ^( m - n ) +2Ae ^>((m|e A:r ( L +- I '-V)) |{A,m})({-A,n}|. (39) 

A m,n 

Noting ((n\e x ^ L+ ~ L ^\n)) = ((n\e~ x ^ L+ ~ L ~^\n)) by the results in section 3 of H F can 
be written as 

H' F = J2 ((n\e x ( L +- L ->\n)} £ e 2iA0W |{A, n}> ({-A, n}|. (40) 
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Here we want to solve the equation i(d / 'dt)^> = H F ^ completely, however it is very 
hard (see for example [S], [U], [211, Therefore let us appeal to a perturbation theory. 

For simplicity we set <p = in (jSZj), then we have the well-known formula 

e 2iX0(t) = £ J a (2\g 2 /u E )e iaUE 1 = J (2g 2 /cu E ) + E J a {2\g 2 /u E )^ E \ (41) 

where J a (x) are the Bessel functions. For a further simplicity we set 2g 2 /ujE = T. 
We decompose (J47IJ) as 

F F = + H lF ; 

where 

^ = E(H^ (i+ - L - ) H)^o(r)El{A,^})({-A,n}| (42) 



and 



4 = E(H^ (VL) I«))EE J »W^1{^}>«-V}I' ( 43 ) 



Next let us transform (JHHJ). 

^e^^Se^^^mle^+^-V)) |{A,m»({-A,n}| 



E e'^^'EE J a (Ar)e^*((m|e Aa:(i+ - L - ) |n)) |{A, m}) ({-A, n}|. (44) 



m,n A Q 



Now we define a new basis called Schrodinger cat states 

\{a^ n }) = ^=(|{l,n}> + ^|{-l,n}», * = ±1- 
From these we have easily 

|{A,n}> = WJ) + A|{-1, V„}», A = ±1. 

Now we are in a position to rewrite (jUj) , (JHSJ), (jSJ) in terms of |{A, VVi}) : 

^F = E((^ (L+ - L - ) H)^o(r)EA|{A,^})({A,^}l (45) 



and 



h' 1f = E ((^|e x(L+ - L - V)) E ^ Et E J > r ) 

n a ^0 a 



2 



j{EA|{A,Vn}}({A,Vn}|-^EA|{A,Vn})({-A,Vn}|}, (46) 
A A J 
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and moreover 



1 

2 



^A|{A,^})({A^n}|-^E A l{ A ^m})({-A,Vn}| ■ (47) 
- A A J 



For simplicity in the following we set 



then 



E 



A,n,A 



^A,n,A = ^A((n|e^ L+ - L -V))^o(r), T = ^ (48 
z we 



(iV) f Ae"^L n (/t 2 ) J (r) where k = x 

(K) f A (2^(1 + /t 2 )-^-"F n (fi: 2 : 2tf) J (r) where « = sinh(x) (49) 
(J) f A^r(l - K 2 ) J - n F n (K 2 : 2 J) J (r) where « = sin(x) 

from the results in sectin 3.1 of |14j . 
Then we rewrite ()35|) as 

and appeal to a perturbation method. Namely, we treat H' 0F a unperturbed Hamiltonian 
and the remainig a perturbed one. 

It is easy to solve the equation i(d/dt)^ = (A/2)H' 0F ^/ , which solution is given by 



^0 = EE e "^ A -' Ac n,A|{A,^}), 
n \ 

where {c n ,\} are constant, so we can set an ansatz to solve (JHOJ) as 

*o = EEe- itEi ^c n , A (i)|{A,^}) (51) 

n A 

and determine the coefficients {c nt \(t)} from (|50p. However the (infinite) equations are 
almost impossible to solve. On the other hand we are interested in Quantum Computation, 
so let us restrict the ansatz (|5T|) which is enough for our purpose. 

One Qubit Case 

n 



The ansatz is simple 

*o= E e- itE ^c n , x (t)\{X^ n }), (52) 
Ae{i -1} 

where n is fixed. Substituting this into ()50|) and after some algebras we obtain 

^\((n\e^- L -^\n)) E e «—*£ (e~^^c n , A (t) - ae"^— c n ,_ A (t 

so we have 
<i 

^c„, A (t) = 



-i|A((n|e^+- L -)|n}) E E (e*— <c n , A (t) - ae^" 2S — ) C „,_ A (t) 

Now we set a resonance condition : for some a^OeZ 

acuB - 2E A;n _ 1 = (-a)u; E - 2E A>riil = 0, (53) 

because E^ n _\ = — -Ea.^a- Then the remaining terms might be neglected (a kind of 
rotating wave approximation), so we have 

By the way, from the fact J- a (x) = (— l) a J a (x) = J a (—x) we have 

5>j_ a 0rr) = 5>j a (-<rr) = E(-^RW = -E^K)- 

cr cr cr cr 

Therefore the equations above become 



d I Cn,l \ _ I *f | [ C n,l 



(54) 



with 



^ = A(He^- L -)|n))E^^^ = A(He- (L+ - L - ) |n))J Q (r) 1 ( o 1} . (55) 



cr 
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We note that this 1Z is just the Rabi (flopping) frequency. The solution is given by 



/ 



c n ,i(t) 



cos(yt) ism(^t) 



\ C n-l(t) 



isin(yf) cos(yt) 




(56) 



Two Qubit Case 



We identify the two-qubit space with two excited states, 0,0, 
For m < n the ansatz is 

*0= E e~ UEA ' m ' X CmAt)\Mm}) + E e- UEA ' n - X CnAt)\{\,^n}). (57) 

Ae{l,-i} Ae{i,-i} 
Substituting this into ()5()j) and after long algebras we obtain 

ie- itE ^j t c m>x (t) = 

^X(( m \e^- L ^\m)) E e**"* E (e- UEA ^c m , x (t) - ae"^— c m ,_ A (t)) + 

A 

~2"~ ^2 



X( jtn(m-n)j2e ia " Et J2 {{m\e axi - L+ ~ L - ) \n)) (f itE ^c n , x {t) - ae- itEA ^c n ^ x (t 



ie- itE ^j t c n , x (t) 
A \{{n\e^+- 1J -'\n)) 



X((n\e^- L ^\n)) E e*—' E (e^^c^f) - ae~ UE ^ c n ,_ A (i)) + 

-Ae lin( "- m) ^ e ^E ((nle^^+^-V)) (e^^c^t) - ae- itE ^-*c m ^ x (t) 



2 V - 

Then we have 



ij tCm , x (t) 



^\((m\e^- L ^\m)) E E (<^*W*) - orf**"" } c m ,_ A (t)) + 



fAEE^y £l ((-|e ra(i+ - i - ) l-))x 



3 ii(ao; B +C(m-n)+£; A>miA -£; A> „ iA ) c _ cre it(aiw s +n(m-n)+BA lTOl A--EA,n,-A) c A (t) 1 
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i Jt c n ,x(t) 
A 



^\{{n\^- L -)\n)) £ E {e^CnAt) - <xe 



_,it(atj_ B -2£; Aini _ A ) 



2 
A 



Cn,-A(0) + 



AEE^^<H e ™ (L+ ~ L_) N> x 

2 ~ ~ 2 



; it(ao; £ ;+C(n-m)+_B A , niA -_E A , miA ) c - (Te it ( m B+ n ( n " m )+ £; A 1 » ! A-EA,m -a) c a(0| 

Now we can set four (possible) resonance conditions : 
(1) a resonance condition : for some a G Z 

OW^ + 0(m — n) + #A,m,l - ^A,n,l = 

(-a)^i? + 0(n - m) + #A,n,i - #A,m,i = 0. 



Then by using the rotating wave approximation we obtain the equations 



d 
dt 



( c 1 ^ 

<^m,l 
Cm, — 1 
Cn,l 

V c « -i / 



4 2 



ft 



£-771,1 

Cm.— 1 
Cn,l 
V C n-1 / 



with 



n = ae ^(He^+^-V)) 

= AJ Q (T)^ {((m|e x(L+ - L - } |n)) + (-l) a <(m| e ~ a,(I ' + "~' L ~ ) 



Here we have used the following identity : neglecting A in 1Z 



n ^ |^^K) ((m | e -(X + -L_) 



^))} t = E 



J (ar) 



((n|e" 



-<rx(L+— L_) 



m}) 



E 



jQ( ~ ffr) ((n|e^- L -) 



m )) = ^ J - a ^ r) ((n| e ^ L +- L -V))- 



2 " 1 1 " V 2 

For the explicit value for matrix element ((m|e CT:r ( i+ ~ L ~)|n)) (a = ±1) see 
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The solution is 



^ ( cos(tf ) 



Cm,l (0 
Cm -lit) 

c n ,i(t) 

\ C n -l(0 J V 

(2) a resonance condition : for some a E Z 



-i_|sin(tM) ^ f Cmi(Q) \ 



^sin(tM) cos(t M) 



c m -i(0) 
c„,i(0) 

V C n,-l(0) 



au E + Q(m -n) + £ A , m ,-i - #A,n ,-i = 
-<=^ (-a)^ + fi(n - m) + #A,n,-i - ^A, m ,-i = 0. 

Then by using the rotating wave approximation we obtain the equations 



d_ 
di 



Cm,l 
Cm,— 1 
C n ,l 



( 



1 2 











Cm,— 1 
Crt,l 



with 

The solution is 



U = AJ2 ^Q((m|e ra ( i +- L -)|n)). 



c m ,i(t) 
c m -i(t) 
c n ,i(0 

V c n _i(t) y 



cos(tf) 



V ^sin(i 



cos(tf ) y 



c m ,i(0) 

Cm,-l(0) 
C n ,l(0) 

V cn,-i(o) y 



(3) a resonance condition : for some a E Z 



ojcjb + fi(m - n) + £ A ,m,i - #A,n,-i = 
(-a)cj s + fi(n - m) + #A,n,-i - #A,m,i = 0. 
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Then by using the rotating wave approximation we obtain the equations 



d_ 

dt 



Cm,l 
Cm,— 1 
C n ,l 

V c «^i J 



( 



o 



% K \ ( 

1 2 



Cm,l 
Cm, — 1 
c rt,l 



with 



K = A^ ( r^^((m|e^ i +- L -V)) 



AJ a (r)i {((m\e xiL +~ L -\n)) - (-l) a ({m\e~^ L +- L -^\n}}} 



Here we have used the following identity : neglecting A in 1Z 



»»} t = E^ 



((n\e 



-ax(L+—L- 



E(-^) J " ( €rT \ (n\e° x{L+ - L - ) \m)) = - ]T a J - a< f T) ((n|e^ L +- L -) \m)) 



so 



E 



(7- 



J - a(ffr) ((n|e-( L +- L -> 



m» = -n. 



The solution is 



c m ,i(t) 
c m -i(t) 

C n ,l(0 

V c n _i(o y 



^sin(tM) ^ f Cml(0 ) \ 



C m -l(O) 

c„,i(0) 

V cn-i(o) y 



V^sin(^) cos(t M) y 

(4) a resonance condition : for some a G Z 

acJe + Q(m - ra) + #A,m ,-i - #A,n,i = 
(-a)u; E + Q(n - m) + #A,n,i - #A,m,-i = 0. 
Then by using the rotating wave approximation we obtain the equations 



d_ 

dt 



Cm,l 



\ 



Cm, — 1 
C„,l 



( 







1 2 



1 2 



W C ! ^ 

Cm,— 1 
C„,l 
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with 



n = A j2 a I^H(( m \e^ L+ - L -^\n)). 



(74) 



The solution is 



^ Cm,l(t) ^ 
C n ,l{t) 



cos(t^) 
-^sin(tf) 



-^sin(tM) 
cos(t^) 



1 c mil (0) ^ 
c m ,-i(0) 

Cn,l(0) 
V C n ,_i(0) j 



(75) 



On the ansatz (jBTj) we solved the Schrodinger equation (|5T))) in the strong coupling 
regime (!) under the resonance conditions and rotating wave approximations, and ob- 
tained the unitary transformations of four types which are a generalization of [Hj. They 
will play an important role in not only Quantum Computation but also Quantum Optics 
or Condensed Matter Physics. 

It is very interesting that each Rabi frequency 1Z contains some Bessel functions { J a (T)}. 
Sec [12J for an interesting experiment which the Bessel function J appeared. Sec also [9_. 

By the way, we considered the case of one atom with two-level, so we would like to 
generalize our method to the case of n atoms (with two-level) interacting both the single 
radiation mode and external periodic fields like (n atoms trapped in a cavity with the 
photon interaction) 

|0) |0> |0> 

^ o o • • • o ^ 

ll> II) ll> 

Then the Hamiltonian may be 

~ n ^ n n 

H nL = tul M ® L 3 + ^ ]T cr? ® (L + + L_) + -J2 ® 1l + 92 E cos^t + ® 1 L , 

j=i j=i j=i 

(76) 
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where M = 2 n and (k = 1, 3) is 

0* = h <£> • ■ ■ <8> I2 <8> crfe <£> I2 ® ■ ■ ■ <S> I2 ( j - position). 

See j21] as an another model similar to this (its model assumes the RWA from the starting 
point). 

In the near future we will attempt an attack to this model. We would like to construct 
C-NOT operators for each pair of atoms (this is a very important subject in realistic 
Quantum Computation). 

By the way, according to increase of the number of atoms (we are expecting at least 
n = 100 in the realistic quantum computation) we meet a very severe problem called 
Decoherence, see for example ^Hj and its references. We unfortunately don't know how 
to control this at the present. 

A generalization of the model to N-level system (see for example ^H], [Uj, pH], [20] ) 
is now under consideration and will be published in a separate paper 1 . 

Acknowledgment. The author wishes to thank Marco Frasca for his crucial suggestions. He 
also wishes to thank the graduate students Kyoko Higashida, Ryosuke Kato and Yukako 
Wada for some help. 
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